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Abstract — In this paper, an improved technique of using 
open-ended rectangular waveguide sensor for 
simultaneous non-destructive complex electromagnetic 
parameters determination of compound absorbing 
materials is developed. The technique is based on 
employing a layered material medium, which consists of 
material to be tested and known materials, sandwiched 
between sensor aperture and perfect conductor plate to 
obtain two complex reflection coefficients necessary to 
achieve this purpose. Finite-Difference Time-Domain 
method is adapted to account for the finite flange size of 
the sensor and numerically calculate the sensor reflection 
coefficient under different test conditions. The technique 
can be employed for cases of single layer and multilayer 
medium planar sheet measurement. The parameters of the 
material under test are obtained by fitting both the 
calculated and measured reflection coefficients using the 
iterative optimization technique. The details related to the 
analysis, FDTD modeling and testing procedure of the 
proposed technique is discussed. The simulations and 
experimental results are compared with the published data 
by literatures and companies to validate the proposed 
technique. The technique is promising for potential 
applications such as design and fabrication of sheet or 
layered coating materials and nanocomposite material and 
bio-medium characterization . 

Keywords — FDTD; Non-destructive test, Dielectric and 
magnetic characterization, layered medium, rectangular 
waveguide sensor. 

I. INTRODUCTION 

Radar absorbing materials (RAMs) is essential in a wide 
range of applications such as radar detection, shielding and 
electromagnetic compatibility technology [1-5]. These 
materials are usually designed and fabricated from 
compound materials as a sheet or coated on metallic 
structures in order to reduce unwanted reflections or to 
reduce the radar cross section area of an object. Moreover, 
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these materials are usually used in a specific application to 
provide a certain performance of reflectivity over a given 
frequency range especially at microwave and millimetre 
wave [6], The design of radar absorbing materials is 
constrained by many factors, including thickness, weight, 
and absorbing properties, which are described by both 
complex permittivity e,= e -j e and complex permeability 
p,= ,u - jp over broadband frequency range. These 
properties are strongly influenced by the frequency at the 
microwave frequency range [7-8]. 

For many years, evaluation of electromagnetic properties 
of materials has been known to be a fundamental aspect of 
microwave and millimetre-wave technology especially 
radar absorbing materials [9-11]. It has aroused 
considerable interests in developing measurement 
techniques to characterize these materials quickly, 
accurately, and conveniently. Recently, various techniques 
are proposed and developed to characterize 
electromagnetic properties of radar absorbing materials at 
microwave frequency range with their respective 
advantages and drawbacks among the others is open-ended 
rectangular waveguide sensor [12-13]. This technique is 
widely used to determine both complex permittivity s r and 
complex permeability p r over a wide band of frequency 
range [14] due to its simplicity and openness in its 
structure. Moreover, it requires less sample preparation. 
These merits make the technique suitable for 
characterizing the electromagnetic properties of materials 
under laboratory and field conditions since the sensor is 
placed directly against the material under test [15]. 
Currently, Finite Difference Time Domain (FDTD) 
method is employed to simulate absorbing ability of radar 
absorbing materials to electromagnetic waves [16]. Also, 
with increasing the use of nanocomposite materials, there 
is a demand for characterization of these materials, starting 
from the designing process up to the fabrication stage [17] 
due to its ability to handle complex structures and 
geometries. In this paper, a technique of using rectangular 

Page | 49 


International Journal of Advanced Engineering Research and Science (IJAERS) 
h ttps://dx. doi. org/1 0. 221 61/ijaers.4. 8. 9 


waveguide sensor irradiating into a layered material 
medium in conjunction with FDTD method is proposed to 
simultaneously determine both complex permittivity and 
permeability (s, and p r ) of compound absorbing materials. 
The structure, which consists of material under test backed 
by known low loss material, is sandwiched between sensor 
aperture with finite flange and perfect conductor plate to 
measure two complex reflection coefficients necessary to 
determine s r and p r . FDTD method is employed to 
theoretically calculate the sensor reflection coefficients 
under different physical test conditions. It is also can be 
evolved to simulate the design of the electromagnetic wave 
absorber. The calculated reflection coefficients of the 
sensor are imposed on the measured values to extract c r 
and p, by iterative -optimization technique over an X-band 
of microwave frequency range. To effectively model these 
materials, the thickness variation and frequency 
characteristics of c, and p r of these materials are studied. 
The FDTD simulations and measurement results on 
selected samples of compound absorbing materials are 
presented and some guidelines for this technique of 
measurement are highlighted. 

n. PRINCIPLES AND ANALYSIS 

Figure 1 shows the geometry of the proposed technique 
used to perform e, and p, measurement. As shown, the 
geometry consists of finite flange open-ended rectangular 



Fig.l: Geometry of the proposed technique: 


waveguide sensor in Cartesian coordinates with wide and 
narrow dimensions (a and b) placed directly against 
layered material medium consisting of the material under 
test with unknown s r and p, and known thickness di, 
followed by known £ r , p, and thickness material and the 
whole structure is backed by a perfect conductor to 
increase measurement sensitivity. The backing known 
material used should be low loss material to increase 
measurement accuracy [18], As shown in Fig. 1, the 
reflection coefficient symbolized r o (a, b, s r , p„ / d) is a 
complicated function of the waveguide sensor dimension 
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a , b, test frequency f sample thickness d, and the 
parameters of the material under test s r and p r . The energy 
radiated from the sensor penetrates through the layered 
material medium and reflected back into the aperture. 
Consequently, the reflection coefficient calculated at the 
sampling point carries the desired information of material 
under test. The EM-parameter extraction of material under 
test by the inverse problem process includes finding the 
values of e r and p, which minimize the difference between 
the reflection coefficient theoretically determined (F t h) and 
the measured one (T m ) obtained under different test 
conditions. Therefore, for accurate extraction of s r and p r . 
the theoretical reflection coefficient of the sensor must be 
accurately predicted. Using the approximate solution for 
theoretical reflection coefficient introduces error in the 
extracted values of s, and p r . Hence, the problem geometry 
is formulated using FDTD method to account for the finite 
size of waveguide flange and accurately predict the 
reflection coefficient under different test conditions. 

1. The Proposed Technique 

The main objective of this work is to develop a technique 
for simultaneous s r and p r determination of compound 
absorbing materials during the design phase and 
fabrication process respectively. In principle, to extract 
four parameters, two complex reflection coefficients are 
required. For measurement purposes, in order to determine 
two complex quantities (c r and p, ) with the assumption 
that the material under test is backed by conductor, it is 
necessary to perform a certain experimental procedure to 
measure the needed two independent complex reflection 
coefficients under two different test conditions. Physically, 
the two conditions can be achieved by different ways. 
Using the proposed technique, the following two ways can 
be performed to obtain these two conditions as shown in 
Fig. 2 ( a and b). 


MUR ABC 



(a) (b) 

Fig. 2: Conditions of reflection coefficient measurement 
(a) Single layer test (b)Two layers test 
1 . 1 Single-Layer Test 
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This test is to be performed to obtain one of the two 
needed reflection coefficients. In this test, the material 
under test with unknown e, and p r and known thickness d\ 
(only single layer) is sandwiched between the rectangular 
waveguide sensor and the perfect conductor as shown in 
Fig. 2 (a). The measured reflection coefficient p m i is 
described using (1): 


Pm\ =T 0 (e r ,/i r ,f ,d) (1) 

As can be seen in (1), the reflection coefficient is 
frequency dependent while the frequency is an 
independent variable. It is to be noted that both complex 
permittivity c, and permeability p, of most materials are 
also frequency dependent and they are slowly varying with 
frequency. 

1.2 Multilayer Structure Test 

This test is to be performed to obtain the second reflection 
coefficient. A medium with two layers is used for this 
purpose as shown in Fig. 2 (b). In this test, the sensor is 
placed in close contact with a multilayer structure 
consisting of the material under test with unknown e r and 
p, followed by known material and the whole structure is 
backed by the perfect conductor. The measured reflection 
coefficient p 2 m is described using (2): 

P m Z= T o( £ r’ £ r \Mr’Pr ’>/ ^ ,d ') (2) 

Where p m i and p m 2 are the reflection coefficients measured 
at different test conditions, £ r , s r , p r and p r are relative 
permittivity and relative permeability of the tested sample 
and known material, d\ and di are thicknesses of material 
under test and known material respectively. It is to be 
noted that using medium with two layers provides enough 
information to extract e r and p r of material under test 
(using both single layer and two layers tests) while using 
only single layer test provides the required information if 
frequency-varying method or thickness -varying method is 
employed. The proposed technique can be also applied to 
perform multiparameter measurement (complex 
permittivity, complex permeability, and thickness) using a 
medium with N different layers. The rectangular 
waveguide sensor used in the measurement has practically 
a flange with finite dimension which makes it possible to 
test high loss sheet samples. For finite flange rectangular 
waveguide sensor, by using the FDTD method, both 
electric and magnetic fields distributions are calculated at 
different physical test conditions. Consequently, the 
reflection coefficient Fo of the sensor can be determined 
according to the theory of transmission-line. 


2. Numerical Modeling using FDTD Method 

The Finite-Difference Time-Domain method has been 
initially proposed by Yee [19] to numerically solve the 
Maxwell equations and developed by Taflove [20]. It has 
been proven to be the most popular numerical technique 
successfully applied for the solution of electromagnetic 
problems with complex geometries [21], In this paper, 
FDTD method is employed to calculate the sensor aperture 
reflection coefficient since it is quite difficult, in this case, 
to apply the analytical approach. Figure 2 (a and b) shows 
the analytical model of the problem in 3D Cartesian 
coordinates considering only the dominant TEio mode of 
the field with sensor dimensions of a and b respectively. 
The boundary space of the problem is divided into two 
regions. The first region is the interior of the waveguide 
(region 1) while the second region is the multilayer 
structure backed by perfect conductor. Each region within 
the geometry is assumed to be linear, isotropic and 
homogeneous in nature and characterized by both £ and 
p. The set of Maxwell’s curls equations are given using (3) 
and (4): 


8t £ 

oH I „ 

= VxE 

8t /u 


(3) 


(4) 


Using 3D Cartesian coordinate system (x,y,z), eqs. (3) and 
(4) are expressed in a system of six coupled Maxwell’s 
differential equation. For x-coordinate of both electric and 
magnetic fields' components, the differential equations are 
described using (5) and (6): 


dH 1 dE dE 

— = ( ~ — ) 

dt ju dy 8z 

dE x 1 dH z dH y 

dt e dy dz 


(5) 

( 6 ) 


Yee has proposed an approach to numerically solve these 
equations by which the problem space is divided into a 
mesh with lattices. A point (i,j,k) in lattice space is 
denoted as: 


(i , j ,k ) = (i Ax , j Ay , k A z ) 
where A is the increment in x, y, z directions. According to 
Yee procedure, both electric field E and magnetic field H 
are distributed to be interleaved with respect to cell whose 
origin is located at (i, j, k). The components of both 
electric and magnetic fields are updated in discrete time- 
steps of At interval. Mathematical theorems for the FDTD 
formulation, concerning issues such as accuracy, 
convergence, computational complexity and stability is 
provided in [20]. Following Yee procedure, the x-direction 
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components of both electric and magnetic fields are 
described using (7) and (8): 



1 



x{i,j,k) 


p Ay 


J7 n — F n 

+ 1) u y(ij,k K 

Az 


(7) 


A 1 H ..-H 

= r + _ { — j ^ — i 

£ Ay 




A z 


( 8 ) 


Following the same procedure, the other field components 
can be obtained and will not repeat here. It is to be noted 
that in this formulation the EM parameter s and p of 
material under test can be constant, or they can be changed 
with frequency. This allows modeling of material 
properties for lossless and lossy over required broadband 
analysis. The total space of the problem is discreetized into 
cubic cells with fine space dimensions to increase the 
accuracy of calculations. To limit the domain of 
computational space, the simulation of open problems is 
usually carried out by placing absorbing boundary 
conditions (ABCs) in the terminating planes of the grid. In 
this paper first-order Mur ABC [22] is used for the 
following reasons: 

- It can be applied to the boundary with a dielectric 
discontinuities 

- As compared to other absorbing boundary conditions, It 
is computationally much easier to implement. 

A forward -moving TEio wave is lunched at the excitation 
plane A- A’. Maxwell’s equations solving algorithm 
involves a continuous sampling of electric and magnetic 
fields in a finite region at equidistance points in a spatial 
lattice. The electromagnetic field component distribution 
(both electric field and magnetic field) is observed within 
waveguide and multilayer structure until the steady state 
condition is obtained. The sampling point (B-B) is chosen 
to be located far away from aperture to avoid higher order 
modes. The complex reflection coefficient T 0 of the sensor 
is calculated by applying transmission-line theory using 

(9): 

r a =K\e- J *- = y\y° (9) 


where Ya is the aperture admittance and Y 0 is the 
characteristic admittances of sensor aperture. 

III. SIMULATIONS AND EXPERIMENTAL 
RESULTS 

1. FDTD Scheme Validation 

To verify and validate the FDTD formulation of the 
problem using the proposed technique, the complex 
reflection coefficients of X-band rectangular waveguide 
sensor irradiating into layered material medium is 
calculated using 3D FDTD code developed for this 


purpose. In this work, we passed the developed code by 
comparing the results of numerical simulations of 
reflection coefficients with those obtained experimentally 
using standard WR-90 rectangular waveguide with flange 
length chosen to be 50 mm [23], The comparison is made 
for two cases; the case when only the material under test as 
shown in Fig. 2 (a) is considered and for the case when a 
medium of two layers is considered as shown in Fig. 2 (b). 
In the first case, the sample used is compound absorbing 
material with s r =14.95 - j.0. 3318 and of p, = 1.35 - j 
1.647 and thickness of 2.8 mm while for the second case, 
the same sample under test is used followed by the second 
layer with known s r and p r and thickness and the whole 
structure is backed by a perfect conductor. Teflon as low 
loss material is chosen [18] to be known material with 
thickness of 8 mm. Simulations and measurements are 
performed at a frequency of 10 GHz using a HP-851 OB 
automatic network analyzer and the results are shown in 
table 1. For the two considered cases, the FDTD 
simulations and experimental results are in good 
agreement validating the computational tool. 

Table 1. Comparison of the sensor reflection coefficient 
obtained using FDTD modeling versus experimental 

results for two cases considered. Results are obtained at 

f= 10 GHz. 


Case 

Method 

in 

(Magnitude) 

T (Phase) 
(Phase in 
Degree) 
dedeDegree) 

Single layer 

FDTD 

0.4391 

-139.38 

Experiment 

0.4424 

-138.17 

Two layers 

FDTD 

0.4062 

-140.28 

Experiment 

0.4132 

-137.98 


2. EM-Parameter Determination 

The analysis based on FDTD modeling in the previous 
sections is to be verified experimentally in order to assess 
the validity of the proposed technique for simultaneously 
determination s r and p r of compound absorbing materials. 
Several experiments are conducted for this purpose over 
the given microwave frequency range. Table 2 shows a 
comparison between the measured results of £ r and p r of a 
compound absorber sample and corresponding previously 
published data in [24] where it is seen that they are fairly 
consistent. The FDTD calculations and measurement are 
performed at a frequency of 10 GHz with sample thickness 
of 2.08 mm. The extraction process of both s r and p, are 
performed iteratively using numerical optimization 
technique. 
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Table. 2: Comparison between the measured results of 
Srand p r of a compound absorbing material and reference 
data Results are obtained at f = 10 GHz. 


Method 

£r=£ -j £ 

pr=p'-jp' 

s' 

e” 

F' 

p" 

Proposed 

Technique 

15.21 

0.37 

1.73 

0.87 

Reference 

Data[24] 

15.26 

0.68 

1.71 

0.84 


As stated in previous sections, the sensor reflection 
coefficient is sensitive to more than one variables with 
operating frequency as independent variable while both c, 
and p, are frequency dependent ( i.e s ,-(/) and g,(/)j. Also, 
the sample thickness is another limiting factor in the 
measurement process. Another test is performed to study 
the influences of these factors on s r and p, measurement. 
The variations of real parts and imaginary parts of both e, 
and p r are separately investigated with frequency and 
sample thickness in this study and the results are illustrated 
in Fig. 3(a and b). Figure 3-a shows the real parts (s and 
p)variations of the s r and p, while the variations of the 
imaginary parts (tan5 E and tand M ) of e, and p r are shown in 
Fig. 3-b. The measurements are performed on another 
compound absorbing material with two thicknesses having 
a ratio of 2. The first sample thickness is chosen to be 2.44 
mm and the second one is of 4. 88 mm thickness. Twelve 
frequency points are selected for the measurement to cover 
the X-band frequency range. It is clear from the two 
figures that good agreement is obtained between the 
measured values (both real parts and imaginary parts for 
sample of 2.44 mm thickness. For sample with 4.88 mm 
thickness, a large discrepancy is observed in the results of 
measurement for both permittivity and permeability. The 
results obtained showed that the sample thickness has an 
influence on measurement accuracy and sensitivity. It is 
quite difficult to obtain accuracy in measurement of 
thicker sample high loss materials such compound 
materials. This is due to. that the reflection decreases with 
increasing the thickness of absorbing material. This is the 
reason why thickness measurement accuracy becomes 
poorer for thicker lossy materials. Hence, backing the 
layered material medium by perfect conductor enhances 
the measurement sensitivity and accuracy. 



8 9 10 11 12 13 

Frequency(GHz) 

(b) 

Fig.2: The measured permittivity and permeability of 
compound absorbing material with tw’o thicknesses 
compared with reference data, (a) Real parts (b) Loss 
tangent. 

IV. PARAMETERS £ r , p r EXTRACTION 

For open-ended rectangular waveguide loaded with 
material under test, the reflection coefficient (or its 
equivalent aperture admittance Y a ) is theoretically 
formulated as a function of sample s r and p r parameters, 
sample thickness d and operating frequency /. As with 
most material characterization methods, there exists no 
closed form, which relates both s r and p, to the theoretical 
reflection coefficient (F thy ). Hence, in order to 
simultaneously extract s r and p, , inverse problem 
technique must be employed numerically for this purpose 
using an iterative optimization technique. By using the To 
value, both obtaining from theoretical analysis and 
experimental under two different test conditions, the 
complex parameters e, and p r of the tested sample can be 
extracted. In this regard, the extracted values of e, and p r 
should be those values taken by which the difference 
between the calculated and the measured values of the 
sensor reflection coefficient is minimized by inverse 
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problem. Consequently, an iterative solution using 
optimization technique is sought to compute e r and p, for a 
given value of reflection coefficient r o to optimize the 
objective function tp using (10). 

^ i (r (/.k qo) 

i . _ i m i c i 

where f is the operating frequency of the measurement 
data set measured at N frequency points. It is clear from 

(10) that by using the values of the measured and the 
FDTD calculated reflection coefficients, the objective 
function (tp) can be optimized iteratively. Minimizing 
process usually starts with initial guess values for e, and p r . 
These iterative algorithms are guaranteed in most cases to 
produce a correct convergence, but this process consumes 
much time. The technique discussed in this paper is easy to 
use. It needs only reflection coefficient measurement, then 
the remaining task is to process the measurement data in a 
convenient way. 

V. CONCLUSION 

An improved technique for simultaneous non-destructive 
EM-parameter of sheet compound absorbing material 
determination using rectangular waveguide sensor 
irradiating into layered medium has been numerically 
investigated and assessed. The FDTD simulations and 
experiments results of both c, and p, variations (real and 
imaginary parts) with the thickness of the sample over an 
X-band of the microwave frequency range have shown that 
the sample under test thickness is strongly influencing s r 
and p, measurement accuracy. The obtained results 
suggest that using the proposed technique for s r and p, 
measurement is suitable for high loss materials with a 
thickness of several millimeters. For this reason, the 
layered structure is backed by a perfect conductor to 
improve measurement accuracy and sensitivity. The 
proposed technique can be applied to test a single layer 
and multilayer material under test where it is quite difficult 
to use the traditional methods such as thickness-varying 
method and frequency-varying method for this purpose. 
The measurement geometry described in this paper is 
limited to two layer structure, which yield enough 
information to determine simultaneously electric and 
magnetic parameters. Since the sensor’s reflection 
coefficient is sensitive to e r and p r and sample thickness, 
the technique can be used to perform simultaneous 
multiparameter measurement (both complex permittivity 
and permeability and thickness). 
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